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We study the magnetic structure of tlie ground state of an itinerant Fermi system of spin- 1/2 
particles with magnetic dipole-dipole interactions. We show that, quite generally, the spin state of 
particles depend on its momentum, i.e., spin and orbital degrees of freedom are entangled and taken 
separately are not "good" quantum numbers. Specifically, we consider a uniform system with non- 
zero magnetization at zero temperature. Assuming the magnetization is along z-axis, the quantum 
spin states are k-dependent linear combinations of eigenstates of the Pauli matrix. This leads 
to novel spin structures in momentum space and to the fact that the Fermi surfaces for "up" and 
"down" spins are not well defined. The system still has a cylindrical axis of symmetry along the 
magnetization axis. We also show that the self energy has a universal structure which we determine 
based on the symmetries of the dipolar interaction and we explicitly calculated it in the Hartree- 
Fock approximation. We show that the bare magnetic moment of particles is renormalized due to 
particle-particle interactions and we give order of magnitude estimates of this renormalization effect. 
We estimate that the above mentioned dipolar effects are small but we discuss possible scenarios 
where this physics may be realized in future experiments. 
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I. INTRODUCTION 

Interest in ultracold atoms with dipole-dipole interac- 
tions arise due to the long range and anisotropic nature of 
dipolar forces. If the dipolar Fermi gas is fully polarized 
the anisotropic form of the interactions leads to a distor- 
tion of the Fermi surface (FS) (Ref. 1, 2 and Appendix 
A) and to a mixing of the collective excitations due to the 
structure of the Fermi liquid parameters.^ Although this 
system can be described in terms of the standard Lan- 
dau theory of the Fermi liquid'^ important changes are 
needed to account for the effects of the dipolar forces. ^'^"^ 
Recent theoretical work has revealed that gases of dipo- 
lar Fermi atomic systems can exhibit interesting and un- 
conventional properties, including novel ordered quan- 
tum liquid crystal states, such as biaxial nematic and 
ferronematic phases^'^. Quantum liquid crystal phases 
have been conjectured to play an important role in the 
physics of strongly correlated systems®. Recent exper- 
imental progress in trapping strongly magnetic dipolar 
Fermi gases in magnetic traps opens the way to study 
these interesting physical systems in the laboratory.^ 

Dipolar Fermi systems where spin degrees of freedom 
are allowed to play a dynamical role"^'^"'^^ are much less 
studied. Magnetic dipole-dipole interactions, unlike elec- 
tric dipole-dipole interactions, conserve total angular mo- 
mentum, spin -f orbital (J = L -|- S),^^, and hence mix 
the spin and orbital degrees of freedom. Interactions in 
dense nuclear matter with non-central forces are simi- 
lar. In fact, many of the above effects have been studied 
theoretically in the nuclear physics literature. 

In this work we consider, in perturbation theory, the 
effects of magnetic dipolar forces on a spin- 1/2 Fermi sys- 
tem with a population imbalance of particles originally 
in eigenstates of the cr^ Pauli matrix, corresponding to 
a partially polarized ferronematic state. ^ To first order 



in the interaction, he local spin quantization axis is not 
along z, leading to an effective spin-orbit type interac- 
tion for the quasiparticles, whose self-energy develops a 
spin structure in momentum space. In an earlier pub- 
lication we showed that the qualitative physics of the 
ferronematic phase can be captured by the structure of 
the Fermi surfaces. The results we present here show that 
the nature of the ordered state at partial polarization has 
a complex, unconventional (i.e., a spin triplet particle- 
hole condensate anisotropic in orbital space), magnetic 
order associated with the dynamically generated spin- 
orbit type interactions. Such a system can be simulated 
experimentally with two hyperfine levels of ^^'^Dy in an 
optical trap at zero external magnetic field. ^ 

We define the dipolar interaction parameter as A = 
rip?' I If, where /i is the magnetic moment of the par- 
ticles, n = fc|n/(37r^) and kp is defined by the average 
chemical potential jl = (/ii +/i2)/2 = ep = h'^kp/2m. 
Our calculations are to first order in A, but the physical 
principles wc describe are quite general. We do not con- 
sider here important questions such as pairing effects, 
collapse instabilities^'*, temperature or optical trap ef- 
fects. The physics we describe does not rely on trap pa- 
rameters; such effects could be taken into account using 
Thomas- Fermi functionals.^^ Collapse instabilities, that 
we will not consider here, require a finite dipolar coupling 
strength even for long range forces^ whereas the phenom- 
ena we study occurs for infinitesimal dipolar couplings. 

Before presenting the details we summarize the physi- 
cal picture that emerges as follows. A fixed difference of 
particles in the eigenstates of ct^ creates a net magneti- 
zation of the system which then acts, as far as static pro- 
cesses are concerned, as a uniform external magnetic field 
along z. In the absence of magnetic field the dipole inter- 
action conserves the spin -I- orbital angular momentum 
(S + L). Crucially, in the presence of a net magnetiza- 
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tion the colliding particles still conserve the z-projection 
of the total angular momentum, = + along the 
magnetization. Recall that the magnetic dipole interac- 
tion is an L = 2 object under rotations of the spatial 
coordinates^^ and that, unlike contact, Coulomb or ex- 
change interactions, the dipolar forces can flip the spin 
of one of the two interacting particles. To conserve an- 
gular momentum the missing momentum is taken from 
the orbital degrees of freedom. In a many-body system 
this means that any given particle will "feel" an average 
magnetic fleld due to the interaction with other particles. 
In other words, the mean field is a vector not a scalar. 
This molecular magnetic field is capable of flipping the 
spin of the spin- 1/2 and hence it must point away from 
the the 2;-axis. Moreover since Jz must be conserved the 
change in spin angular momentum due to a spin flip pro- 
cess must be accompanied by a corresponding change in 
orbital (momentum space) angular momentum. In this 
way we arrive at the conclusion that the distribution in 
momentum space of the particles must be anisotropic and 
that the spin state of the particle will depend on its ve- 
locity! 

Since the; local spin quantization axis is not parallel to 
the polarization direction z, the k-dependent spin state 
of the particles becomes a linear combinations of eigen- 
states of cr^, i.e., the spins tilt. Interestingly, we find that 
neither the tilting angle nor the momentum distribution 
of the quasiparticles are isotropic in momentum space, 
and vary with direction as ^ |F2±i(k)| and l2o(k) respec- 
tively because of the mixing of orbital and spin degrees 
of freedom. This can be described in perturbation the- 
ory as arising from two processes: a) |t) — >■ |t) (no spin 
flip), with ALz = 0; since we consider particles with no 
internal structure, the quasiparticlc self-energy acquires 
a correction of the form F2o(k), b) |t) — > \i) (spin flip) 
with ALz = TI7 the quasipiarticle self-energy gets a cor- 
rection of the form |y2.±i(k)|. Therefore the self-energy 
now has an off-diagonal momentum-dependent compo- 
nent. These changes in the structure of the self-energy 
lead to qualitative modifications to the quasiparticlc dis- 
tribution functions and of the Fermi siirface. In par- 
ticular it is no longer possible to define separate Fermi 
surfaces for up and down fermions independently. 

This work is organized as follows. In section II we in- 
vestigate the structure of the single particles states and 
show how the orbital and spin degrees of freedom become 
entangled due to the presence of non-diagonal elements 
in the self energy of particles. In section III we show that 
particle interactions lead to a renormalization of the bare 
magnetic moment. In section IV we compute the oc- 
cupation number in momentum space noting that this 
quantity is a 2 x 2 tensor in spin space. In section V we 
provide possible scenarios where dipolar effects may be 
observed. Finally the appendices provide the details of 
the self energy computations for the fully polarized and 
partially polarized case (Appendix C and D), an intu- 
itive explanation of the FS distortion (Appendix A) and 
details of the Fourier transform of the dipolar interaction 



(Appendix B). 

II. STRUCTURE OF THE QUASIPARTICLE 
SELF-ENERGY 

The excitations of a gas of fermions at finite density 
are quasiparticles with the same quantum numbers as the 
fermions. However their dynamics is strongly affected by 
their mutual interactions. These physical effects are cap- 
tured by the quasiparticle self-energy. In the presence 
of dipolc-dipole interactions the form of the self energy 
is determined by symmetry arguments. We thus require 
that the quasiparticle self energy be invariant under si- 
multaneous rotations of the net magnetic field and the 
quasiparticle momentum, which leads to an expression 
of the form 

I]^^(k,u;) = MK^)S^j + ^^^^(3kik, - dij) (2.1) 

where Sq and arc functions of the wave vector and 
frequency, i = x,y, z is the direction of the effective mag- 
netic fleld and j = x,y,z is the direction of the spin 
polarization. The form of this expression is valid to all 
orders in perturbation theory. 

Without loss of generality, we consider a microscopic 
Hamiltonian with a population imbalance of particles in 
eigenstates of a^, i.e. we will assume a state with a fixed 
total polarization. The Hamiltonian is 

k,a k 

+ ^T.S\mj{ci)S^{-ci) (2.2) 
q 

where Ck = are the singe particle energies, nk,t '^nd 
nk,i are the occupation numbers of the single particle 
states with momentum k with both spin projections, and 
g is the strength of the repulsive short-range s-wave in- 
teraction. The last term is the dipolar interaction. In 
Fourier space it is given by (see Appendix B) 

V,,{q) = (4V/3)(3q.qj - 6,,) (2.3) 

where q is the unit vector in the direction of q, /x = 
is the bare magnetic moment of particles and 7 the 
gyromagnetic ratio, e.g., for electron spins, 7 = 2/is. 
Repeated indices are to be summed over. 'S'*(q) = 
X^k ^k+q a'^afj'^KIi IS the Spin density in momentum 
space, cr* are Pauli matrices and a, /3 are spin indexes 
that label the eigenstates of a^. 

In what follows we will assume that either the system is 
prepared in a state of finite polarization or that g is large 
enough for the ground state to be a ferronematic with 
finite polarization (as discussed in Ref.[7]), and compute 
the fermion self-energy perturbatively in the dipolar in- 
teraction. To first order, the self energy is given by (see 
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Fig.l) 



(Pk' 1 ^ 



X gO^(k',*w„O^.j(k-k')aL5<0 (2-4) 
Explicitly (see appendix D) 

h^^2 = = -/(fc)3k,(k, + zky). (2.5) 

Being a 2 x 2 matrix we can expand ?iI]Q,^(k) in Pauli 
matrices, KEap = S^'(k)tT^^, where we found Y,^^{k) — 

— /(A;)(3kik2 — 5i^z)- The angles of k are defined in Fig. 
3. Comparing with the general expression Eqn. (2.4), we 
conclude that in the Hartree-Fock approximation 
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St - -2/(fc). 



(2.6) 
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FIG. 2: (color online) characteristic form of /(fc) (blue-solid) 
and its derivative (red-dashed) for two values of /ir and p., 
Eqn. (2.7). Main panel: kpi ~ 1.1, kF2 = 0.9, inset: kpi = 
1.8, kF2 ~ 0.3 in units of kp- 
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FIG. 1: Feynman diagrams for Hartree (left) and Fock (right) 
terms of the self energy 



We assumed a fixed population imbalance of particles 
in the two eigenstates of , by introducing two chemical 
potentials fXi and ^2- G°^(k,'ia;„) = Sa/s/iii^n - Cka) 
is the free temperature Green function, where ft^ka — 

~ A*Q is the bare single particle dispersion relation and 
^k = f?y? I'lm. In Eq.(2.6) /(fc) is a smooth monotonic 
function of the magnitude of the wave vector: 



TT 



dr 



j2{kr)jn(k'r) 



x[npiei,-^il)-nF{ei,-^i2)] (2.7) 



where npix) = [exp (/3a;) -I- 1]^^ is the Fermi function. 
The self energy vanishes for an equal number of particles 
in each eigenstate of a^. This means that dipole torques 
exist only when there is a net magnetization in the sys- 
tem, as we would expect. 

The diagonal matrix elements of the self energy, pro- 
portional to ±(3cos^ — 1), give anisotropic renormal- 
ization to the single particle energies. It was show that 
for partial polarization the majority-component FS is 



elongated while the minority-component FS compression 
along the polarization axis.^ The novel feature of Eqn. 
2.5 is its non-diagonal matrix elements. They are not 
zero because the dipole interaction can flip one of the two 
interacting particles. This means that the mean held act- 
ing on a given particle's spin can flip it. In turn this imply 
that such an effective molecular vector field is not along z- 
axis. The structure of the terms: Ei2(k) ~ y2,-i(k), and 
E2i(k) ~ l2,i(k) of the non-diagonal matrix elements 
can be understood as a consequence of angular momen- 
tum conservation. Since the magnetic dipole-dipole in- 
teraction conserves Jz = Lz + Sz in the presence of a net 
polarization along z-axis, a change in the spin angular 
momentum ASz = ±1 in accompanied by a correspond- 
ing change in the orbital angular momentum ALz = Tl- 
The particle energy, (ek)a^, have a tensorial struc- 
ture and is given by the equation h^uaSais + ^^^^(k) = 
(ek)Q^ — Expanding in Pauli matrices we obtain 



(ek)Q^ 



■ ^k '^ai3i 



(2., 



we see that there is an the effective magnetic field, hk, 
given by ej^ = fih]^, that acts on the particle's magnetic 
moment. This k-dependent field has a contribution from 
the external field long z (i.e., fj,r) and a molecular con- 
tribution due to interactions with other particles. This 
field is necessarily not along z. We explicitly write the 
components e]^ in Fig. 3 and we plotted hk in Fig. 4. 



III. QUASIPARTICLE MAGNETIC MOMENT 

As we remarked above the problem of a system with a 
fixed population imbalance is equivalent to have a static 
magnetic field in the system. Hence it is of interest to 
know what is the magnetic response of the system. In 
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FIG. 3: (color online) Effective magnetic field, (red arrow), at 
point k in momentum space. The local spin quantization axis 
is tilted with respect to z because the dipole-dipole interac- 
tion induce non-diagonal matrix elements, Ei2(k), in the self 
energy, Eqn. (2.5). 



a small external field we define the effective magnetic 
moment of the particle as 



/iij(k) 



^4 



5eL 



(3.1) 



Where we defined the external field along z by ^lhl^^ = 
/ir . The effective magnetic moment of the particles is not 
the bare magnetic moment ^. In fact, when dipole-dipole 
interactions are present we can guess the form of the 
magnetic moment carried by the particle by symmetry 
arguments. We require that the magnetic moment be 
invariant under simultaneous rotations of the magnetic 
field and the momentum of the particle, i.e.. 



AT(k,w) 



(3k,k, -,5„-) (3.2) 



where we find jlo — M and jiT — 2/i(i9/(fc)/9/ir). f{k) is 
given by Eqn. 2.7. For the values used in Fig. 4 we find 
P't/P'O ^ 0-3 at the outer FS. For small dipolar couphngs, 
where our microscopic calculations are accurate, the oc- 
cupied states in momentum space lie outside the singular- 
ity, see Fig. 4. The line singularity is a ring in the kx,ky- 
plane with radius, k^, given by /ir — f{ko). This wave 
vector, fco ^ {iir/ fJ-'^Y^^ , is perpendicular to the polar- 
ization axis, assumed long z, and naturally introduces a 
length scale, 1 /ko into the system. In the regime were our 
calculation is valid the length scale is much sorter than 
the inter-particle distance, 1/ko < l/m.eLx{kpi,kF2) and 
hence it is related to microscopic effects not considered 
by our long-distance microscopic dipole-dipole Hamilto- 
nian. For strong dipolar couplings, our calculation sug- 
gests that the radius of the ring penetrates the FS's. In 




ky 



FIG. 4: (color online) 2D view of points in k-space where 
Ek± — fi{red and blue). FS's of free system: = fi ± 
/^r (dashed lines). Superposed are stream lines of magnetic 
moment field carried by particles in momentum space, Eqn. 
3.2. This field resembles the magnetic field produced by an 
electric current circulating counterclockwise in a ring placed 
in the horizontal plane. Values for these plots are kpi = 1.1 
and kF2 = 0.9, in units of kp, and A = 0.3 << 1. 



this case, particles would occupy spin states with non- 
trivial topological structure. At a critical coupling the 
system may eventually acquire a spatial modulation. 
The strongly coupled regime, A — > oo, can be accessed 
by increasing the dipole interaction or by reducing the 
kinetic energy. The latter possibility can be experimen- 
tally engineer in the laboratory by placing the atoms in 
an optical lattice. 



IV. ONE-BODY DENSITY MATRIX 

In this section we compute the one-body density ma- 
trix in Hartree-Fock approximation. The distribution 
function of particles, (7ik)a/3, is now a 2 x 2 matrix in 
spin space. Dyson equation, {G^^)ai3 — {i^^n — £,'ka)Sai3 ~ 
I]c(^(k, uj), with Hartree-Fock self energy corrections pro- 
vide a conserving approximation^^ to the Green function. 
Defining fi.^k± = Ek± — an explicit calculation gives 



522(k, iWn) 

fJi2(k, iw„) 

f/2l(k, iuJn) 



1 



Uk^ke 



1 



iw„ - ^k- 
1 



(4.1) 



where Ek± = ± [(?iEii(k) - ^,)2 + |ftEi2(k)|2] and 
U]^ — cos(6k/2), i'k = sin(6k/2) and < 6k < tt is de- 
fined in Fig. 3. is the probability of up-spin propaga- 
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tion and is the probability of down-spin propagation. 
Equating Q = QqI + Qia^ + G2'^^ + Gsc^^ and collecting 
like terms we can write in compact form 



{iuJn - ^k+)(«Wn - Ck-) 



(4.2) 



where ^kM = (^k+ + Ck-)/2 = (e[^ - p,)/h and ^kr = 
(^k+ — Ck-)/2. hk is a unit vector along the effective 
magnetic field parametrized by the polar angles 0k,</>k, 
see Fig. 3. 1 is the 2x2 identity matrix. In our 
case the one-body density matrix is given by (nk)a/3 = 
J(c?a;/27r)ni?(a;)/3Q^(k, oj), where pa/B is a 2 x 2 tensor of 
spectral functions and np(u!) = [exp (/?tj) + 1]^^. The 
spectral functions are given by 

-^(k,-")|..„=.+.o+ (4-3) 

In Hartree-Fock, we obtain two simple delta functions 
corresponding to the poles of the Green functions. It is 
useful to expand the density matrix in Pauli matrices, 



(4.4) 



where rik = (l/2)tr[(nk)] is the average occupation 
of state with momentum k, and (<S'^) = ^"-k = 
(fi./2)tr[cr*(rik)] is the i-th spin density component. (5^) 
has a clear physical interpretation namely, it is a vector 
field with polar angles defined in Fig. 3 and plotted in 
Fig. 4, which gives the k-dependent quantum mechanical 
average of the spin state of the particle. Explicitly, 



n-k 

(-Sk) 
(SI) 



-nF(Sk-)] 

h Uk^k cos (/)k [np (£^k+ ) - np {Ei^- )] 
/lUk^kSin^k [nF{Ek+) - TipiEi^^)] 



(Si) = -iul-vl)[npiEi,+ )-nFiE^^)]. 

Note that in these expressions we define np{x) — 
[exp {f3x — fl) + l]^^. Since ^ and ^ 0, there 
is a non-zero average probability that the spin is point- 
ing along the xy-plane in the region between the blue and 
red curves in Fig. 4. One can check that 



0. 



(4.5) 



i.e., there is no overall spin magnetization of the gas in 
the xy-plane. The system has cylindrical symmetry along 
the magnetization axis, see Fig. 4. 

The net spin magnetization is along z and is given by 
Sk ("^k) ■ "Th^ average occupation number of state k and 
spin projection along z: rik^ = n]^+n^ and rik^ = n\^—nf^, 
are 



nw,t = ulnp{Ek+) + vlnpiEk-) 
nk,l = Wk"F(-E'k+) + ulnpiEi^-)- 





FIG. 5: (color online) density plot of the occupation proba- 
bility in momentum space of ^-projection of spin, nkt (a) and 
nki (b), (Eqn. (4.6)). Values for these plots are the same as 
those of Fig. 4. 



Fig. 5 is a density plot of the occupation probabilities of 
z-component of the spin in momentum space. Note that 
the average occupation of spin states is not uniform in 
momentum space as it would in a non-interacting Fermi 
gas in which case there would be two concentric spheres 
uniformly occupied with the "up" or "down" eigenstates 
of . Note also that dipole-dipole interactions cause the 
spins to tilt predominantly along the diagonals, see also 
Fig. 4, in such a way that measurement of quantum spin 
state will yield "down" with non-vanishing probability 
even in regions in momentum space originally occupied 
only by spins pointing up. 

Eqn. (4.6) implies that the ground state wave function 
which evolves adiabatically from the non-interacting sys- 
tem is of the form 



I* 



n ( 



keOk 



|0) 



(4.6) 



where Uk is the k-dependent amplitude for the spin to 
be the eigenstate "up" of cr^ and Wk the k-dependent 
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FIG. 6: (color online) momentum distribution along a line 
with polar angle 9^ = 60 deg (see Fig. 3) for rik^ (red-solid) 
rik^ (blue-dashed) . 



amplitude for the spin to be the eigenstate "down" . In 
this work we computed these coefficients to first order 
in the dipole interaction. As we see the spin states of 
particles are linear combinations of eigenstates of a^. 
The true FS, the one which separates the occupied from 
unoccupied single particle states, is given by the red 
curve in Fig. 4 i.e., by the values of that satisfy 

^2 , ifcv.^n,M2li/2 



E^,+^fl = el^+[{hJ:n(k)-^^rr+\h^M\^^''~^ = o. 

There are sharp discontinuities in momentum space, see 
Fig. 5. In this sense they define elliptically distorted 
FS's such as those found in Ref.[7]. We comment that as 
the polarization of the gas increases decreases and for 
a fully polarized gas we recover the elliptical shape, red 
curve in Fig. 4 found in Refs. [1,2,18]. 



V. DISCUSSION 

We have used perturbation theory to find the mag- 
netic structure of the ground state of an itinerant, mag- 
netized, spin-Y2 Fermi gas with dipole forces. Our cal- 
culations are the first step in the analogous perturbative 
calculations of Gell-Mann and Briickner of the energy of 
electron gas and hence are of general theoretical interest 
in themselves. However, we have ultracold atomic gases 
with strong magnetic moments, such as -^^'^Dy (lO^s) 
as a specific candidate^^. The dipole interaction of two 
magnetic dipoles in a Fermi gas with density 10^'^ cm""^ 
is ~ lOnK and may be too small compared with other 
energy scales in the system. We expect that it is an 
experimental challenge so observe directly dipole effects. 
In fact, direct evidence of this unconventional magnetism 
may become possible using a variation of NMR experi- 
ment on ultracold dipolar gases as in the superfluid phase 
just as in '^He-A phase^". We briefly consider the follow- 
ing scenarios of possible realization of dipolar effects in 
Fermi systems. 

The "spin" states could be two hyperfine levels of 
i63j)y g^jj optical trap. In its ground state, "'^^^Dy has 
spin, orbital, nuclear and total electronic angular mo- 
mentum of 5' = 2, L = 6, / = 5/2 and J = 8 respec- 



tively. The total angular momentum of the atom can be 
F = 11/2, ...,21/2. However, the energy gaps are several 
orders of magnitude larger that the typical dipole-dipole 
energy per particle attained with current experimental 
densities and its possible that the spin structures we de- 
scribe here are washed out. However, since we are op- 
erating in an optical trap, it may be possible to use the 
low field limit from the same hyperfine-Zeeman manifold, 
in which case, the splitting between hyperfine-Zeeman 
states is not large. Now, if we operate in the high field 
limit (as we may want Og = 0, the Feshbach resonance in 
Dy) then clearly the effects of dipole interaction are sup- 
pressed. At the time of writing little is known about the 
theory Feshbach resonances in Dy. In a closely related 
problem we assumed that the imbalanced Fermi gas was 
metastable. A possible stabilizing mechanism is provided 
by a Stoner-type of effective interaction^^'^^ which was 
not considered in the present work. Another possibility 
is to choose two hyperfine states for which conservation of 
angular momentum would strongly limit the phase space 
available for particles to decay into. If this is true for Dy 
then there is no need to invoke the Stoner mechanism. 

A second possibility is the use of an optical lattice 
loaded with hyperfine states of ^^'^Dy. This method 
has certain advantages as far as stability is concerned^'^ 
and might also make the dipolar effects much more pro- 
nounced. By working with a small filling fraction to avoid 
a large same-site interaction and umklapp processes it 
may be possible to make the dipole interaction be the 
dominant energy scale 

Another important experimental question is the pos- 
sible presence of external stray magnetic fields. In the 
present work we already assumed a fixed external Zee- 
man field. Hence, the same considerations apply as re- 
gards to the observability of dipolar effects under stray 
magnetic fields. In particular if < Ed, the external 
field will pin the overall axis of symmetry of the sys- 
tem but will not modify significantly the spin structures 
described here. We only considered zero temperature 
effects. We known that temperature effects decrease the 
FS distortions^'^^"^^ . In general we expect the spin struc- 
tures we found to be very sensitive to temperature. For 
example, for current attainable densities(10^'^CTO~'^) the 
dipolar energy per particle is Ed ^ lOnK. Boundary ef- 
fects are also known to be very important for long range 
forces such as dipolar. If a degenerate dipolar Fermi gas 
is realized one will obtain images that resemble Fig. 5 in 
time-of-flight experiments. 
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Appendix A: The distortion of the Fermi surface of 
a one-component dipolar Fermi gas 

One of the most interesting theoretical predictions is 
a FS which is eUiptical in shape with symmetry along 
the magnetization of the system. Such predictions have 
been made for a dipolar one component Fermi gases^'^'^* 
and recently for two component Fermi gases'". The fact 
that the ground state FS have elliptical shape in fully po- 
larized systems can be understood in several ways a) by 
calculation of the self energy corrections to the single par- 
ticle energy levels^, see appendix C for details, b) by not- 
ing that if the density is kept constant it acts as a source 
field which entangles itself with the quadrupolar mode^. 
c) by noting that the symmetry of the bare interactions 
in a fully polarized system is invariant under rotations 
about the net magnetization axis. An elliptical FS with 
with same axis of symmetry is consistent with the sym- 
metry is a possible ground state. This argument, how- 
ever, does not tell whether the FS is prolate or oblate^. 
Here we provide perhaps the simples qualitative argu- 
ment to understand the effect. It is based on elementary 
kinematics of a particle in an anisotropic potential and 
the Pauli principle. Consider a gas uniformly polarized 
along 2:-axis. The Pauli principle prevents occupation 
of the same single particle states and hence gives some 
stability to the system.^'' The interaction between two 
dipoles is V{r) = (/i^/r^)(l — Scos^^r), where 9r is the 
angle of the relative position vector with the z-axis. As 
we see, the energy of particles in a region close to the xy- 
plane(6'r ~ n/2) experience, on average, a repulsive inter- 
action, i.e., they move "up the hill" . The total energy per 
particle is h"^ (kxy) /2m -|- n/x^ = E. The mean separa- 
tion between particles is n~^/^ and hence = n, the 
particle density. However, particles in the spatial region 
close to the z-SLxis{9j. ^ 0) experience on average an at- 
traction, i.e., they move "down the hill" and hence the to- 
tal energy per particle is fi^ (k1) /2m — 2n/i^ = E. Com- 
paring this two expression we find for particles with same 
total energy (fc^) > (k^y). Of course (fc^) = {kxy} = 
because of translational invariance. This suggests that 
the kinetic energy and hence the Fermi wave vector 
of particles close to the z-axis is bigger than the ki- 
netic energy of particles close to the xy-plane. Defining 
Kpx = {k'pxv) /2"^ and Kpz = (^Fz) /2?7i from the 
relations shown above we obtain 

Kfz-Kfx ~ 3nd^. (Al) 

As we would expect, the difference in the Fermi energies 
along the z-axis and xt/-plane is proportional to the dipo- 
lar energy per particle. For spin-independent forces in a 
translational invariant system the wave vector of parti- 
cles, k, is a "good" quantum number. The FS, which 
concerns only with this quantum number is decoupled 



from the spin state of the particles. Surprisingly, dipole- 
dipole interaction furnish an example where this is no 
longer true and k is entangled with the spin. 



Appendix B: Fourier transform of the dipolcir 
interaction 

The Fourier transform of the dipole-dipole interaction 
is obtained by Fourier transforming the left hand side of 
the identity 

= ^(<5i, -3fif,), (Bl) 

and using the well known results that the Fourier trans- 
form of l/|x — x' I is 47r/g^ . Here r = x — x', r = |x — x'| 
and f is a unit vector along r. The contact term does 
not contribute for finite separation between dipoles but 
is required because the diagonal elements ( i = j) on 
the left hand side must satisfy the Poisson equation 
-V2(l/|x-x'|) = 47r(5(x-x'). This makes the left hand 
side and right hand side traceless. It is also required 
because the magnetic field is divergence free. One could 
see this property directly from Eqn. (2.3). Note also that 
the dipole interaction is only valid at long distances and 
hence here we only consider the physics of a system of 
particles whose separation is larger that the Bohr radius, 
i.e., the electronic clouds of the atoms do not overlap. 



Appendix C: Self energy of fully polarized dipolar 

gas 

In this section we consider a homogeneous Fermi gas 
interacting with dipolar forces. The results of this section 

apply to electric as well as magnetic dipoles. The calcula- 
tion below was outlined in Ref . [2] . Using the formula for 
the decomposition of a plane wave into a sum of spherical 
harmonics we can show that the Fourier transform of the 
dipolar interaction Vd(r) = (/x^/r'^)(l — 3cos-^ Or), can be 
written as 

= Vlm,Vm'{k,k')Yi*^{V!)Yi,^,{ii) 

lm;l'm' 

(Cl) 
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Where 

Vim,i'm'{k,k') = -2(f{Anf{-i)¥ 



X (^J ^jiikr)ji'{k'r) 

J dnY;^{v)Yi'm'{^)P2m 



(C2) 

The Hartree contribution to the self energy vanishes be- 
cause Vd(q = 0) = 0. The Fock term can be written 
as 



Si/F(k') 



d^k 



V{k-k')nl 



(27r)3 

= ^CTim{k)YUk') 



Im 



(C3) 



Where we defined 



(^l,m{k) = j ^^Vlm,Vm'{k,k')Yiim'{^)nl 

I'm' ^ ' 

(C4) 

In perturbation theory we use the spherical distribution 

in momentum space. We calculate the term Voo,;m 
and substitute into Eqn. (C4) to obtain 

crim = -(>l,20mfl \/ / k dk 

TT V 5 7o 



/■°° dr 

/ -k{kr)j2{k'r) (C5) 
Jo r 



Where k^ is the undistorted(spherical) Fermi wave vec- 
tor. Wc conclude that the Hartree- Fock self energy is 
P2(k), explicitly 



Sff^(k) =a(/e)P2(k) 



(C6) 



where cr(fc) is a smooth slowly varying function of the 
magnitude of the wave vector 



a{k) 



367rA;3 



'iek"p + ^k-'{k''pf -ik{k%~f 



-3(A;'- (4)^)^arctan(P) 
k> . 



(C7) 



This equation is valid for all k with fc< being the least 
of k'p and k. One can check that cr(fc) is continuous and 
twice differcntiable at k^p. Wc conclude that the single 
particle dispersion to first order is given 



= —+a{k)P2{cos9k) 



(C8) 



In particular the Fermi wave vector can be written as 



ki? = nkp + nSkp 



(C9) 



To first order in the dipolar interaction the chemical po- 
tential is unchanged, n = {kp)^ /2m. The first correction 
to the chemical potential is quadratic order in the inter- 
action. The Fermi wave vector is given 



k2 



Solving to first order in Skp we obtain 



(CIO) 



Skp = a^f^L^s cos^ 0^ _ 1) (Cll) 



Which means that the FS is elliptical in the weak cou- 
pling limit for fully polarized dipolar gas. 



Appendix D: Self energy of an imbalanced 
two-component dipolar Fermi gas 



In this section we compute in Hartree-Fock the self 
energy of a polarized spin-i/2 system. Starting from the 
expression 



?lSa/3(k) 



/ 



d^k' 1 ^ 



xGg^(k', ^Wr^')Vij(k - k')<,a^^, (Dl) 



one can show that S21 = ^^'^ ^22 = — Perform- 
ing the spin indexes summation and the Matsubara sums 
we obtain 
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;iSn(k) = - 
;iSi2(k) = - 



3 



/drk' 
^^3cos6lk-k' sinfilk-k/e-"^"-"' [nF{ei, - /ii) - npiei, - ^2)] 



(D2) 



r 



Note that these expressions vanish for a system with 
equal number of particles in each of the eigenstates. 
By expanding the dipolar interaction in spherical har- 
monics and using the unperturbed momentum distribu- 
tion of particles , the integrals can be calculated analyt- 
ically with the result 



;ii]u(k) = -/(fc)(3k2-i) 

fiSi2(k) = -/(fc)3k,(k, -i^) 



(D3) 



where f{k) is a smooth function of the magnitude of the 
wave vector, 

f(k) = ^ Hk'^dk' r -j2{kr)jo{k'r) 

Jo Jo r 

X [npiei, - Ml) - npiei, - M2)] (D4) 

The self energy can be expanded in Pauli matrices as 
S = Sol + Sio-^ -I- S2(t2' + SsCT^ and we obtain 



/1SJ2 — 



(D5) 



where Si(k) = — /(fc)(3kjk^ — Si,^)- Wc know give an 
explicit derivation of Eqn. (D3). Starting from Eqn. 
(D2) and expanding in spherical harmonics 

y2o(k-k') = l!^(3cos2^k-k'-l) 

= EE^2o"^'''"'('^'^')^™(k')^''™'('^) 

Im I'm' 

(D6) 

Substituting back into Eqn. (D2), we obtain the form of 
Ell is 



?lSii(k) = - (71' m' {k)Yi,ra' (k), 



(D7) 



where ai>m'{k), depend only on the magnitude of the 
wave vector, k, 



(7l'm'{k) = ^ 



Im 



v^-''"^\k,k')Y;^ik') 



(27r) 



3 "^20 



X [npiei, - 111) - npiei, - ,,2)] (D8) 

The angular integration is simplified because the un- 
perturbed distribution function is spherically symmetric. 



We obtain 

c^vm'ik) = / k''dk'V°,'''"''{k,k') 

X [npiei, - Ml) - npiei, - M2)] (D9) 

The terms V^^'' (fc, k') can be found by noticing that 
the Fourier transform of V2o(r) = (p^/i''^){^ — 3cos^^r) 
is V2o(q) = (47r/x^/3)(3cos^ 6'q — 1) and by definition one 
can write 

1/20 (fc -k')= J drh-'^^-^'>''V2oir) (DIO) 
Note that V{r) can also be written as 



V2o{r) = -^2J-^Y2o{f) 



(Dll) 



and substituting the spherical wave expansion of a plane 
wave 

e-*- = 4nY,{-iyjt{kr)Yim{m*m{r) (D12) 

Im 

into Eqn. (DIO) we can read the only non-vanishing co- 
efficient 



00;20 
20 



{k, k') = (47r)>2 2 



-Y 



00 



dr 



h{kr)jo{k'r). 
(D13) 

Substituting back into Eq. (D9) we obtain the only non- 
vanishing coefficient 

<T2o(fc) = 2y^^M' j k'^dk' j ^ j2{kr)jo{k'r) 
X [npiei, - Ml) - npiei, - 1^2)] 



-m. 



(D14) 



At zero temperature the unperturbed Fermi surfaces are 
sharply defined by wave vectors kpi = (67r^ni)^/^ and 
kF2 = (677^712)^/^ and f{k) can be analytically computed 
as /(fc) — g{k, kpi) — g{k, kp2) where y) is given by 



9{x,y) = 



7Q-K 



3x°y + 8x^y^ - 3xy 



-^ix'-yrin 



y + x 



y-x 



(D15) 
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Finally, If /ii > /i2 then kpi > kp2 and f{k) > 0. A similar 

calculation yields the quoted result for S12. 

;iSn(k) = -/(/c)(3cos2 0k - 1) = -/(fc)(3k^ - 1) 

(D16) 
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